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We analyze in the framework of pQCD the properties of the single spin asymmetry in heavy flavor 
production by linearly polarized photons. At leading order, the parallel- perpendicular asymmetry 
in azimuthal distributions of both charm and bottom quark is predicted to be about 20% in a wide 
region of initial energy. Our analysis shows that the next-to-leading order corrections practically 
do not affect the Born predictions for the azimuthal asymmetry at energies of the fixed target 
experiments. Both leading and next-to-leading order predictions for the asymmetry are insensitive 
to within few percent to theoretical uncertainties in the QCD input parameters: mg, fiR, fiF, Aqcd 
and in the gluon distribution function. We estimate also nonperturbative contributions to azimuthal 
distributions due to the gluon transverse motion in the target and the final quark fragmentation. 
Our calculations show that nonperturbative corrections to a _B-meson azimuthal asymmetry are 
negligible. We conclude that measurements of the single spin asymmetry would provide a good test 
of pQCD applicability to heavy fiavor production at fixed target energies. 
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I. INTRODUCTION 



In the framework of perturbative QCD, the basic spin-averaged characteristics of heavy flavor hadro-, photo- and 
electroproduction are calculated up to the next-to- leading order (NLO) During the last ten years, these NLO 

results have been widely used for a phenomenological description of available data (for a review see [^). At the same 
time, the key question still remains open: How to test the applicability of QCD at fixed order to the heavy quark 
production? On the one hand, the NLO corrections are large; they increase the leading order (LO) predictions for both 
charm and bottom production cross sections approximately by a factor of 2. For this reason, one could expect that 
the higher order corrections as well as the nonperturbative contributions can be essential in these processes, especially 
for the c-quark case. On the other hand, it is very difficult to compare directly, without additional assumptions, 
the pQCD predictions for spin-averaged cross sections with experimental data because of a high sensivity of the 
theoretical calculations to standard uncertainties in the input QCD parameters. The total uncertainties associated 
with the unknown values of the heavy quark mass, mg, the factorization and renormalization scales, fip and fj,R, 
^QC D and the parton distribution functions are so large that one can only estimate the order of magnitude of the 
pQCD predictions for total cross sections ||^,||. In particular, at the energies of the fixed target experiments the 
theoretical calculations are especially sensitive to the value of the heavy quark mass which plays the role of a cutoff 
parameter for infrared singularities of the theory. For the shapes of the one- and two-particle differential distributions 
the above uncertainties are moderate in comparison with the ones for total cross sections; however they are also 
significant. In fact these uncertainties are of the same order as the contributions of nonperturbative effects (such as 
the primordial transverse motion of incoming partons and the hadronization phenomena) which are usually used for 
a phenomenological description of the charm and bottom differential spectra . 

Since the spin-averaged characteristics of heavy flavor production are not well defined quantitatively in pQCD it is 
of special interest to study those spin-dependent observables which are stable under variations of input parameters 
of the theory. In this paper we analyze the charm and bottom production by linearly polarized photons, namely the 
reactions 
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j]+N^Q{Q)+X. (1.1) 

In particular, we calculate the single spin asymmetry parameter, A{s), which measures the parallel-perpendicular 
asymmetry in the quark azimuthal distribution: 

I da{s,tp = 0)-dcr{s,ip^7r/2) 

dCT(s, = 0) + d(7(.s, = 7r/2) ■ ^ ■ ^ 

Here da{s,Lp) = ^{sj'p) , i s the degree of linear polarization of the incident photon beam, ^/s is the centre of 
mass energy of the process (y_j) and f is the angle between the beam polarization direction and the observed quark 
transverse momentum. The main results of our analysis can be formulated as follows: 



-At fixed target energies, the LO predictions for azimuthal asymmetry (1.2) arc not small and can be tested 
experimentally. For instance, 

A{s = 400 GeV^) |L0^„„« A{s = 400 GeV^) |L0^„„,« 0.18. (1.3) 

-The NLO corrections, both real and virtual, practically do not affect the Born predictions for A{s) at fixed target 
energies. 

-At energies sufficiently above the production threshold, both leading and next-to-leading order predictions for 
A{s) are insensitive (to within few percent) to uncertainties in the QCD parameters: mq, ^jj, /xj?, Aqqd Siiid in the 
gluon distribution function. This implies that theoretical uncertainties in the spin-dependent and spin-averaged cross 



sections (the numerator and denominator of the fraction (1.2), respectively) cancel each other with a good accuracy. 



Our analysis shows that the pt- and x i?-spectra of the single spin asymmetry also depend weekly on the variations 
of QCD parameters. To be specific: Theoretical uncertainties in the pQCD predictions for the azimuthal asymmetry 
differential distribution are always significantly smaller than the ones in calculations of the shape of the corresponding 
spectrum of unpolarized cross section. 

We conclude that the single spin asymmetry is an observable quantitatively well defined and rapidly convergent in 
pQCD. Measurements of asymmetry parameters would provide a good test of the fixed order QCD applicability to 
charm and bottom production. 

Another important question is how the fixed order predictions for asymmetry parameters are affected by non- 
perturbative contributions usually used for the description of unpolarized spectra. We estimate the nonperturbative 
corrections to A{s) due to the transverse motion of partons in the target. For this purpose we use a parametrization of 
the intrinsic transverse momentum distribution proposed in ref. | [lO| (so-called fc^-kick). It is shown that the fc^-kick 
corrections to the 6-quark azimuthal asymmetry A{s) are negligible. Because of the smallness of the c-quark mass, 
the fc^-kick corrections to A{s) in the charm case are larger; they are of order of 20%. 

Modeling nonperturbative effects with the help of the kx-idck and the Peterson fragmentation function we analyze 
also the px- and xp -behavior of azimuthal asymmetry in D- and B-meson photoproduction. It is shown that for the 
B-meson both the fc^-kick and the fragmentation contributions are small almost in the whole region of kinematical 
variables px and xp. For the Z?- meson, the fci^-kick corrections are essential at low values oi px, Pt ^ ™ci and rapidly 
vanish at px > rric. 

The paper is organized as follows. In Sect. II we analyze the properties of heavy quark azimuthal distributions at 
leading order. We also give the physical explanation of the fact that the pQCD predictions for A{s) are approximately 
the same at LO and at NLO. The details of our calculations of radiative corrections are too long to be presented in this 
paper; they will be reported separately in a forthcoming publication ||l^. In Sect. Ill we discuss the nonperturbative 
contributions caused by both fc^-kick and Peterson fragmentation mechanisms. A comparison of the QCD predictions 
with the Regge model ones is also given. 

II. SINGLE SPIN ASYMMETRY IN PQCD 

A. Leading order predictions 

At leading order, 0{aemCts): the only partonic subprocess which is responsible for heavy quark photoproduction is 
the two-body photon-gluon fusion: 
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-fik^) + gikg) ^ QiPQ) + Qip-). 
The cross section corresponding to the Born diagrams (see Fig.l) is: 



dxd(p 



{s,i,Lp,fi%) = C 



em 



2(1 -/32)(/32 _i:2) 



l-x 



Z.2 



(l-i2)2 



(1 +'P-yCos2(p) 



(2.1) 



(2.2) 



where 7^^ is the degree of the photon beam polarization; is the angle between the observed quark transverse 
momentum, pq.t, and the beam polarization direction. In (2.2) C is the color factor, C — Tp =Tr(r°T'')/(A''2 — 1) = 
1/2, and eg is the quark charge in units of electromagnetic coupling constant. We use the following definition of 
partonic kinematical variables: 



s — (k^ + fcg) ; 
u = {kg -pqf ; 



t= {k^- pq) ; 

i — rrin 
x = l + 2 — ^ 



(2.3) 



In the Born approximation, the invariant cross section of the single inclusive hadronic process, 

-f{k-,) + N{kN) ^ Q{pq) + X{px), 

can be written in the form 



EQd?a 



(s, t, u,ip) = j dz6 {s + i + u- 2mQ) g (z, ^^"^ {s, i, (fi, 



(2.4) 



(2.5) 



Here g (z, y^p) describes the gluon density in a nucleon N evaluated at a factorization scale ^p- The hadronic variables 
are related to the partonic ones as follows: 



s — (fc^ + ki^y — s/z] 
t = (/c^ -Pqf = i. 



u=(kN-pQ) ='u/z-m|(l/z-l); 



(2.6) 



In this paper we will discuss photoproduction processes only at fixed target energies. For this reason we do not 
take into account the contribution of so-called hadronic or resolved component of the photon. This contribution is 
small at the energies under consideration Q. 

To analyze the azimuthal distributions in photoproduction, it is convenient to use the parallel-perpendicular asym- 
metry parameters. Apart from the quantity A[s) defined by (1.2), we will consider also the parameters ipr) 
A' {xp), 



a^{pI) 



1 d2o-(s, = 0) - d^ (T(s,py, = 7r/2) 
V-y d^a{s,p\,ip = Q) + A^a{s,p\,'^ = 1^/2)'' 



(2.7) 



1 d^a{s,xp,ip = 0)-d^(7[s,xp,ip^TTl2) 
^^■^^ V^d^cj{s,xp,v = Q) + A^<j{s,xp,^ = T:/2y ^ ' 

which describe the dependence of azimuthal asymmetry on the transverse momentum, p^ = pq rp, and on the Feynman 
variable, xp = pi/pi^max, of observed particle, respectively. The quantities d?(T{s,p\, ip) and d^a{s,xp, ip) are the cross 
section (^.S] ) integrated over xp and over p^, respectively. The quantity dcr(s,<y9) in ( |l.2| ) corresponds to the cross 
section ( |2.5D integrated over both xp and p\. 

Unless otherwise stated, the CTEQ3M [|l^ parametrization of the gluon distribution function is used. The default 
values of the charm and bottom mass are rric = 1.5 GeV and rrib = 4.75 GeV; A4 = 300 MeV and A5 = 200 MeV. 
The default values of the factorization scale /.ip chosen for the A{s) asymmetry calculation are ^p \ charm= 2mc for 
the case of charm production and ^p \ Bottom = nib for the bottom case P,[lO|. Calculating the pp- and xp -dependent 
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asymmetries (p^) and A'^ {xp), we use /ii? Icharm~ S-^/m^ +py for charm and lBottom~ \A"^+Pt ^^o'^ bottom. 
For the renormahzation scale, iij^, we use f^n = jip. 



Let us discuss the pQCD predictions for the asymmetry parameters defined by (O), (2/7) and (p.S|). Our calcula- 
tions of A{s) at LO for the c- and fe-quark are presented by solid lines in Fig. 2 and Fig. 3, respectively. One can see 
that at energies sufficiently above the production threshold the single spin asymmetry A{E^) depends weekly on 
= [s — m\i)/2mN- 

The most interesting feature of LO predictions for A{E~i) is that they are practically insensitive to uncertainties 
in QCD parameters. In particular, changes of the charm quark mass in the interval 1.2 < rric < 1.8 GeV affect the 
quantity ^(^^7) by less than 6% at energies 40 < i?^ < 1000 GeV. Remember that analogous changes of rric lead to 
variations of total cross sections from a factor of 10 at i?^ = 40 GeV to a factor of 3 at iS^ = 1 TeV. The extreme 
choices m?, = 4.5 and mi, = 5 GeV lead to 3% variations of the parameter A{E^) in the case of bottom production 
at energies 250 < E^ < 1000 GeV. The total cross sections in this case vary from a factor of 3 a,t E^ = 250 GeV to a 
factor of 1.5 at E^ — 1 TeV. The changes of ^(£'7) are less than 3% for choices of /ip in the range ^nii, < fxp < 2mb. 
For the total cross sections, such changes of ^p lead to a factor of 2.7 at E^ = 250 GeV and of 1.7 at £'7 = 1 
TeV. We have verified also that all the CTEQ3 versions of the gluon distribution function ||l^ as well as the CMKT 
parametrization | |l3| lead to asymmetry predictions which coincide with each other with accuracy better than 1.5% 0. 

In the Born approximation scaling holds: with a good accuracy the quantity A{s) is a function of the variable 77, 
rj — ArriQ/s, so that 



^(■s) lBottom~ ^( — I*) Icharm ■ (2-9) 

mi 



-b 



The scaling behavior of ^(77) (i.e. its independence of Aqcd/itT'q) is demonstrated in Fig.4g. 

In Fig.5/Fig.6 we show by a solid line the LO predictions for the p|.-/pT-dependcnt parameter A^ (^pp)/AP (pt) 
describing the azimuthal asymmetry in charm/bottom production, pr-distributions grow rapidly with px up to 
Pt ~ ruQ and fall slowly as pr becomes larger than the heavy quark mass. The dependence of py-spectra on 
variations of mg is shown in Fig. 7 and Fig. 8. One can see that, contrary to A{s), the form of the p|n-spectrum is 
sensitive to the value of charm quark mass. Note however that variations of the A^ (p^) shape due to changes of rric 
are significantly smaller than the corresponding ones for the unpolarized distribution ^-{stPt) We have also 

verified that a change in the other input parameters has practically no effect on the computed values of A^ (px) for 
both charm and bottom. 

In the Born approximation the pr-distribution of azimithal asymmetry is practically a function of two variables: 
rj = Attiq/s and ^ = p^/rnq. The function A^ (^) at different values of 77 is shown in Fig. 9. 

We have calculated also the dependence of single spin asymmetry on xp. The LO predictions for A^ {xp) defined 



by ( |2.8D in the case of c- and 5-quark production are presented by solid lines in Fig. 10 and Fig. 11, respectively. 



Xi?-distributions take their maximal values approximatelly at xp « 0.5 and fall to zero as xp ^1. 



B. Radiative corrections 



Let us briefly discuss the NLO corrections to the single spin asymmetry A{s). It is well known that, at order 
C'(c^em 0^5)1 the main heavy quark photoproduction mechanism is the real gluon emission in the photon-gluon fusion: 

7(fc7) + gikg) ^ Q{pq) + Q{p-) + g{pg). (2.10) 



'^Note however that the fixed order predictions for azimuthal asymmetry are insensitive to uncertainties in QCD parameters 
only at energies sufficiently above the production threshold. For instance, at E~, < 10 GeV, changes of the charm quark mass 
in the inte rval 1.2 < rric <1.8 GeV lead to 100% variations of the QCD predictions for A{s). 

^Eq.(L3) does not contradict to eq.(^.9|). As one can see from Fig. 4, the asymmetry parameter A{ri) is not a monotonic 
function of variable 77. So, eq. ( p^ reflects only the fact that A{ri = 0.25) A{ri = 0.025). 
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According to the analysis HQ, practically whole contribution to the heavy quark production in photon-hadron 
reactions at fixed target energies, Ej 1 TeV, originetes from the so-called initial-state gluon bremsstrahlung (ISGB) 
mechanism which is due to the Feynman diagrams with the i-channel gluon exchange (see Fig. 12). Since the gluon 
distribution function increases very steeply at small z, g (z) cx 1/z, the order-a^ photon-hadron cross section is 
determined at i?^ ^ 1 TeV by the threshold, s ^ 4:mQ, behavior of the •yg cross section. Near the threshold, a 
large logarithmic enhancement of the <-channel gluon exchange diagrams contribution takes place in the coUinear, 
Pg,T 0, and soft, pg 0, limits Effectively, the ISGB contribution is proportional to the Born jg differential 

cross section: 

J -ISGB ^ j^LO 

^^{s,t)^^Kis)^{s,t), s^4ml (2.11) 

where K (s) is an enhancement factor. 

Since the azimuthal angle ip is the same for both 7g and QQ centre of mass systems in the collinear and soft 
limits, it seems natural to expect that the eq.( ^.ll| ) can be generalized to the spin-dependent case substituting the 



spin-averaged cross sections in (2.11) by the (^-dependent ones. Indeed, the threshold enhancement of the ISGB 
mechanism is due to the gluon propagator pole in the diagrams in Fig. 12 which is a common factor for both spin- 
dependent and spin- independent amplitudes. 

It is well known that the shapes of differential cross sections of heavy quark production in photon-hadron reactions 
are not sensitive to radiative corrections 1^,^. For instance, the normalized quantity 

where aint{s) = / dp|,^^(s,pf.) and p\ = Pq^t^ is practically the same at LO and at NLO. This fact means that non- 
small contributions of the enhancement factor K (s) (NLO corrections) to the pr-dependent and to the pr-integrated 



photon-hadron cross sections, da^^^ {s^p^) « ^dzg{z)d(j^^'^^{zs,p'^), cancel each other in the ratio ( 2.12 ) with a 
good accuracy. One can assume that the same situation takes place also for the single spin asymmetry which is a 
ratio of the (/^-dependent cross section to the (^-averaged one. 

Our calculations show that it is really the case. Radiative corrections, both real and virtual, to the photon- 
gluon fusion mechanism practically do not affect the Born predictions for the single spin asymmetry in heavy flavor 
photoproduction at fixed target energies. Moreover the azimuthal asymmetry is independent (to within few percent) 
of the theoretical uncertainties in the QCD input parameters {niQ, ^r, iip and Aqcd) at NLO too. The details of 
our NLO analysis are too long to be presented here and will be given in a separate publication 

As to the photon-quark fusion, jq — > QQq, the dominant production mechanism in this reaction is the so-called 
flavor exitation (FE), also arising from the diagrams with the i— channel gluon exchange in Fig. 12. However the 
contribution of the FE mechanism may be essential at superhigh energies only For instance, the contribution 
of the 75 reactions to the unpolarized bottom photoproduction makes only 5-10% from the contribution of the jg 
processes at ^ 1 TeV. It is evident that an account of the photon-quark reactions cannot affect significantly the 
predictions of the photon-gluon fusion mechanism at fixed target energies in the polarized case too. 

III. NONPERTURBATIVE CONTRIBUTIONS 

Let us discuss how the pQCD predictions for single spin asymmetry are affected by nonperturbative contributions 
due to the intrinsic transverse motion of the gluon and the fragmentation of produced heavy quark. Because of the 
c-quark low mass, these contributions are especially important in the case of charmed particle production. In our 
analysis, we use the MNR model |lo| parametrization of the gluon transverse momentum distribution, 

kg — ZftN + kg^T- (3.1) 

According to [lO[ |, the primordial transverse momentum, fcg,T, has a random Gaussian distribution: 



1 d^A^ 1 ^ kj. 



(3.2) 
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where = k^rp. It is evident that the inclusion of this effect results in a dilution of azimuthal asymmetry. In [p|,p^, 
the parametrization (|]^) (so-called fcT-kick) together with the Peterson fragmentation function |l^] , 

D{y) = ^2 ^ (3.3) 

y[l-l/y~s/{l~-y)f 

(where is fixed by the condition j^AyD{y) — 1), have been used to describe the single inclusive spectra and the QQ 
correlations in photo- and hadroproduction at NLO. It was shown that the available data on charm photoproduction 
allow to choose for the averaged intrinsic transverse momcmtum, (/cy), any value between 0.5 and 2 GeV^. 

Our calculations of the parameter A{s) at LO with the fc^-kick contributions are presented in Fig. 2 and Fig. 3 by 
dashed ((fc^) — 0.5 GeV^) and dotted {{k"^) ~ 1 GeV^) curves. One can see that in the c -quark case the /cT-kick 
reduces the value of A{s) approximately by 15% at (/c^) = 0.5 GeV^ and by 20% at (fc|,) = 1 GeV^. The fcy -kick 
corrections to the bottom asymmetry are systematically smaller; they do not exceed 5% for both (fc^) = 0.5 GeV^ 
and {kl) = 1 GeV2. 

Nonperturbative contributions to the pT-dependent asymmetry parameters {pr) ^-nd A^ (pt) are shown in Fig. 5 
and Fig. 6. It is seen that the fc-r-kick corrections decrease rapidly with the increase of the heavy quark transverse 
mass. For the i3-meson, they are negligibly small in the whole region of px- In the case of a Z? -meson production, 
the fc^-kick corrections are essential only in the region of low pT ^ rric. 

Calculating the hadronization effect contributions we use for the parameter e that characterizes the Peterson 
fragmentation function the values Ed = 0.06 for a ZJ-meson and Eb = 0.006 for a B-meson ||l^. Strictly speaking, 
according to the factorization theorems, the application of the fragmentation function formalism can only be justified 
in the region of high px- It is seen from Fig. 5 and Fig. 6 that for px > mq the account of the fragmentation function 
( |3.3[ ) leads to a reduction of the p^-spectra. Fragmentation corrections to A^ (Pr) ^''^^ '^^ order of 15% in the case of 
a charmed meson production; they are less than 5% for a _B-meson case. 

For completeness, in Fig. 10 and Fig. 11 we presented nonperturbative contributions to the xj^-distributions of single 
spin asymmetry. One can see that the fcy-kick corrections to A^{xf) in the bottom case are small in the whole 
region oi xp- As expected, the Xi?-distribution of the c-quark azimuthal asymmetry is more sensitive to the /c^-kick 
corrections. 



So, we can conclude that nonperturbative corrections to the &-quark asymmetry parameters ( |l.2| ) and (2.7) due to 
the fc^-kick and Peterson fragmentation effects practically do not affect predictions of the underlying perturbative 
mechanism: photon-gluon fusion. 

To illustrate how strongly the azimuthal distributions depend on the basic subprocess dynamics let us consider 
the mechanisms of the photon fusion with nonvector primordials. The latter can be nonperturbative, color or white 
objects (say, reggeons at pj^ « 0). In the Born approximation, the partonic cross sections corresponding to the 
reactions 



7 + O ^ Q + Q (3.4) 



have the following form: 



dxdip s 



(l-i:2)2 



{l+'P-yC0S2ip) 



(3.5) 



dxdip s 1 — ' 

where the cross sections a'-^^^ and a'^^^ describe the photon fusion with a scalar (O — S) and a pseudoscalar 



(O = P) object, respectively. One can see from ( |2.2| ) and (3^), (3^) that only the gluon contribution leads to a 
positive azimuthal asymmetry. In the case of a scalar primordial, the asymmetry parameters are predicted to be 
negative. The photon fusion with a pseudoscalar object (say, Tr^-meson) leads to an -independent cross section. 
Note that the expressions ( t3.5|) , (|3.6| ) do not correspond to any concrete models; they are presented only as illustrative 
exercises demonstrating how the azimuthal distribution reflects the nature of the underlying subprocess. 

Completing the Section let us compare the QCD predictions for a £)-meson single spin asymmetry with the Regge 
pole model ones. In principle, at arbitrary values of variables xp and pr, the spin dependent part of a Z?- meson 
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photoproduction cross section can be of two types. The first, | [e, nojprl^ = Pt ^'^^^ corresponds to the contribution 
of the t-channel exchanges with natural spin-parity while the second, Kepr)!^ = Pt cos^ (p, describes the contribution 
of the t-channel singularities of unnatural spin-parity. In the above expressions, e is the polarization vector of the 
photon, uq is a unit vector in the direction of the photon momentum and ip is the azimuthal angle of the Z?-meson 
momentum about ng: tamp = {[e, nolpT) / {epx) ■ Really, the Regge pole model is only applicable in the narrow 
region xp ~ 1 and <C s (so-called triple reggeon exchange limit). In this region, the cross section of the inclusive 
reaction is determined by the contribution of the leading (rightmost in the j-plane) Regge trajectory. The diagram, 
which describes the pseudoscalar meson photoproduction in the triple reggeon exchange limit, is sketched in Fig. 13. 
Since the relevant leading Regge trajectory, ao*, is of natural spin-parity, the spin structure of the corresponding 
reggeon-particle vertex, jaD*D, is written as: 

\P-tao-'D{PT)\'^ OC |[e,no]pT|^ = ^Pt (1 - COS 2^) . (3.7) 

It is seen from ( [3.7| ) that, in contrast to QCD, the Regge model predicts in the region ~ 1 a large negative 
azimuthal asymmetry such that A^{xp) |f^°sgc^ _i ^.^ _^ i 



IV. CONCLUSION 



In the investigation of mechanisms which are responsible for heavy flavor production at fixed target energies, the 
reactions (1.1) with a linearly polarized photon beam are of special interest. As it is shown in the present paper, the 
QCD LO predictions for the single spin asymmetry in both charm and bottom production are non-small and can be 
tested experimentally; in a wide region of initial energy the parameter A{s) defined by ( |l.2D is of order of 20%. Our 
analysis shows that the NLO corrections practically do not affect the Born predictions for A{s) at fixed target energies. 
Unlike the unpolarized cross section, the single spin asymmetry is an observable quantitatively well defined in QCD. 
In particular, at both LO and NLO, A{s) is independent to within few percent of the theoretical uncertainties in toq, 
fiR, fip, ^QCD and in the gluon distribution function. Our calculations show that the pT- and a^^^-distributions of the 
azimuthal asymmetry in bottom production are practically insensitive to nonperturbative contributions due to the 
primordial transverse motion and Peterson fragmentation effects. 

So, the single spin asymmetry in heavy flavor production by linearly polarized photons is an observable quantitatively 
well defined, rapidly convergent and insensitive to nonperturbative contributions. Measurements of the azimuthal 
asymmetry in bottom production would be a good test of the conventional parton model based on pQCD. 

Due to the c-quark low mass, nonperturbative contributions to the charm production can be essential. As it is shown 
in our paper, the pQCD and Regge approaches lead to strongly different predictions for the single spin asymmetry 
in the region of low pT and xp ~ ^- Data on the pT- and Xi?-distributions of the azimuthal asymmetry in I?- meson 
production could make it possible to discriminate between these mechanisms. 

Concerning the experimental aspects, the linearly polarized high energy photon beams can be generated using the 



Compton back-scattering of the laser light off the lepton beams (see, for instance, [16 1^ and references therein). Ac- 
cording to the above references, this method promises to provide the beams of real photons with a definite polarization 
and high monochromaticity. 
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FIG. 1. Photon-gluon fusion diagrams at leading order. 
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FIG. 2. Single spin asymmetry, A{E^,), in the c-quark production as a function of beam energy E-y 
QCD LO predictions with and without the inclusion of fer-kick effect. 
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FIG. 3. The same as in Fig.2, but for the 6-quark case. 
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FIG. 4. Scaling behavior of the asymmetry parameter A{ri) as a function oi rj = AniQ/s at QCD LO. 
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FIG. 5. pr-distribution of azimuthal asymmetry, A''{p^), in a D-meson production; the QCD LO predictions with and 
without the inclusion of the fcr-kick and Peterson fragmentation effects. 




FIG. 6. The same as in Fig. 5, but for the pr-distribution of azimuthal asymmetry, A''{pt), in a B-meson production. 
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FIG. 7. Dependence of the py-spectrum, A^{p^), on the c-quark mass, rric, at QCD LO. 
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FIG. 8. Dependence of the pT-spectrum, A''{pt), on the 6-quark mass, mj, at QCD LO. 
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FIG. 9. Scaling behavior of the asymmetry parameter j4^(C) as a function of ^, ^ = p^jm^Q-, at different values of r\ at QCD 
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FIG. 10. KF-distribution of the single spin asymmetry, J^(xf), in a D-meson production; the QCD LO predictions with 
and without the inclusion of the fcT-kick and Peterson fragmentation effects. 
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FIG. 11. The same as in Fig. 10, but for the case of a B-meson production. 
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FIG. 12. f-channel gluon exchange diagrams in the photon-gluon fusion. 
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FIG. 13. Triple reggeon diagram for the process — > DX. 
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